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It is found that the average baryon phase-space density 
obtained from the ratio of deuteron to proton yields is nearly 
constant over centrality in Au+Au collisions at the AGS. The 
finding ofi^ers an explanation for the puzzling centrality inde- 
pendence of the ratio of charged kaon total yields. The corre- 
lation between the charged kaon ratio and the average baryon 
phase-space density is studied for central heavy-ion collisions 
of various systems over a wide range of beam energy. It is 
found that the charged kaon ratio and the average baryon 
phase-space density both increase with decreasing beam en- 
ergy, and are strongly correlated. Such study may provide a 
new approach to search for medium efi'ects on the kaon mass. 

PACS number(s): 25.75.-q, 25.75.Dw 



Strangeness production has been extensively stud- 
ied in heavy-ion collisions because enhanced strangeness 
production may signal the formation of Quark-Gluon 
Plasma. It has been observed that kaon production per 
participant nucleon increases with centrality in Si-|-Al, 
Si-|-Au, and Au-|-Au collisions at the Alternating Gra- 
dient Synchrotron (AGS), reaching a factor of 3-4 en- 
hancement in central collisions with respect to p+p in- 
teractions 01 . However, the ratio of charged kaon total 
yields [K^ /K^) varies little with the collision central- 
ity [0 . Similar results have been also observed in Ni-|-Ni 
collisions at the Schwer lonen Synchrotron (SIS) ||] 
and Pb-f Pb collisions at the Super Proton Synchrotron 
(SPS) 1^. The nearly constant /K~ is puzzling be- 
cause and K~ are thought to be produced by differ- 
ent mechanisms: K~^s are produced by pair production 
together with a i4r+, while if+'s can be produced, in ad- 
dition, by associate production together with a hyperon. 
These different production mechanisms lead to different 
rapidity distributions which are observed, namely, the 
K'^ rapidity distribution is broader than K~'s [Q. One 
naively expects that in heavy-ion collisions the relative 
contribution of associate over pair production increases 
with centrality, because the associate production thresh- 
old for kaons is lower than the pair production thresh- 
old, and there are more particle re-interactions (at lower 
than the full energy) in central than peripheral collisions. 
Therefore, one expects an increasing K'^ /K~ with cen- 
trality. 

On the other hand, the constituent quark model [||J|] 
has been successfully applied to describe particle ratios 
in heavy- ion collisions |qj7[]. In the constituent quark 



model, = us and — us. Hence, /K' depends 
on the baryon (baryon— antibaryon) phase-space density 
established in the collision zone at chemical freeze-out 
(when particles cease to interact inelastically hence parti- 
cle abundances and ratios are fixed) In this picture, 
the observed constant /K~ implies a constant baryon 
phase-space density over the collision centrality. The ex- 
traction of the baryon phase-space density at chemical 
freeze-out heavily relies on chemical equilibrium models; 
such model studies have been so far limited to central col- 
lisions. On the other hand, the baryon phase-space den- 
sity at final kinetic freeze-out (when particles cease to 
interact both inelastically and elastically) can be read- 
ily extracted from deuteron coalescence data. This is 
because deuterons, due to its small binding energy, are 
mostly formed by a proton and a neutron overlapping in 
phase-space, and cannot subsequently survive collisions 
with other particles. In this Letter, we study the kinetic 
freeze-out average baryon phase-space density ((/)s'") 
a function of centrality, especially for Au+Au collisions 
at the AGS. We show that the value in Au+Au col- 
lisions is nearly constant over centrality. The result may 
indicate that the average baryon phase-space densities at 
kinetic and chemical freeze-out ((/)§'') are connected in a 
way that is independent of centrality. We further study 
the correlation between /K~ and in central 
heavy- ion collisions and demonstrate that / K~ in- 
creases with from SPS, to AGS, to SIS (decreasing 
beam energy). 

However, there is a physics origin that may result in an 
opposite behavior, a decreasing /K^ with increasing 
{D'b^j and that is mass modification in nuclear medium. 
It is predicted that the effective mass is lower in nu- 
clear medium than in free space, and the effective 
mass slightly higher ||l^,|ll[|. A decreasing K~ (and/or 
an increasing K^) effective mass would yield a reduced 
/K^ ratio. In fact, the KaoS Collaboration has stud- 
ied /K~ in Ni+Ni collisions at the SIS and inferred a 
significant drop in the K~ effective mass [1). On the the- 
ory side, however, different conclusions can be reached 
depending on model assumptions p^ , p^ . In this Letter, 
we take a different approach, i.e., studying K'^ / K~ as 
a function of (/)§'". If the mass modification is large, 
one may observe a spectacular phenomenon: K'^ /K~ 
increases and then decreases with (/)§'". Such phe- 
nomenon is not observed, however, we argue that this 
new approach may provide a more direct way to search 
for kaon medium modifications. 
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The Letter is organized as follows. First we present the 
rapidity distributions of the average nucleon phase-space 
density at the AGS for collisions of various centralities. 
Then we extract (/)§'" for each centrality and show that 
the value of is nearly constant. Finally we demon- 

strate that K'^ I is strongly correlated with (/)§'" 
using the central collision data and discuss implications 
of the results on medium effects to the kaon mass. 

In the coalescence model, the neutron phase-space den- 
sity is related to the ratio of deuteron to proton yields 
(d/p) [p^~p^. Following Ref. |l6| , but not assuming iden- 
tical neutron and proton yields, we obtain the spatial av- 
erage nucleon phase-space density at kinetic freeze-out 
as 



1 + p/n 1 idN/dy)d 
2 3{dN/dy)p' 



(1) 



where p/n = {dN / dy)p / {dN / dy)n, and {dN/dy)p.n,d are 
the rapidity densities of proton, neutron, and deuteron, 
respectively. The first term arises from isospin asymme- 
try and will be referred to as 



ai = (l + p/n)/2. 



(2) 



The value of p/n depends on the combination of the pro- 
jectile and target species, the collision centrality and nu- 
cleon rapidity. The value of p/n may vary little with nu- 
cleon transverse momentum because of the isospin sym- 
metry of strong interaction. 

The deuteron and proton data from Si-f Al and Si-|-Au 



at 14.6 AGeV/c |17| and Au+Au at 11.6 AGeV/c |1J] 
are used to extract {f{y))^™- We assume p/n = 1 
for the nearly isospin symmetric Si-|-Al collisions. For 
Si-f Au and Au-|-Au, the values of p/n are estimated using 
the Relativistic Quantum Molecular Dynamics (rqmd) 
model jl^] . For the rapidity ranges covered by the data, 
the model indicates 0.8 < p/n < 1, hence 0.9 < aj < 1 
for both systems. 

Figure ^ shows the obtained {fiy))^™ as a function 
of the proton rapidity, (/(y))^'" varies little with rapid- 
ity in central Si-f Al and Au-f Au collisions. In periph- 
eral Si-|~A1 and Au+Au and both central and periph- 
eral Si+Au collisions, {f{y))^"^ is larger at target rapid- 
ity than mid-rapidity. We note that at target rapidity 
deuterons may be produced by mechanisms other than 
coalescence (e.g. target fragmentation) |18 , resulting in 
an overestimate of (/(y))]^'" by Eq. @). The magnitude 
of this possible effect is not investigated in the present 
work. 

The average baryon phase-space density is the relevant 
quantity. Hyperons are the most significant contributors 
to the baryons besides nucleons. The abundance of an- 
tibaryons is negligible at AGS energies, and is small at 
SPS energies. Generally, the average baryon phase-space 
density is obtained from 

" = (1 - N/N)[l + {Y- Y)/{N TV)], (3) 



where (/)]^™ is the rapidity averaged nucleon phase-space 
density (see below), and N, N, Y, and Y are the nucleon, 
antinucleon, hyperon, and antihyperon total yields, re- 
spectively. It has been assumed that the phase-space 
distributions of hyperons and antibaryons are the same 
as of nucleons . This assumption does not introduce 
significant errors on as the contributions of hyper- 

ons and antibaryons are generally small (see Table |). 
We denote 



a^v = 1 - N/N « 1 -p/p 



(4) 



where N /N has been approximated by the antiproton to 
proton ratio (p/p), and 



Hence, 



ay = 1 + iY -Y)/{N - N). 



{f)B = {f)N ■ O^N ■ ay. 



(5) 



(6) 



The conservations of global baryon number and 
strangeness give, respectively. 



N-N = iVpart -{Y-Y) 



(7) 



whereiVpart is the total number of participant nucleons, 
and 



Y -Y K K - K = {I + uk) ■ {K+ - K-) 



(8) 



where uk = [K^ - K^)I{K+ - R-). Therefore, Eq. (|) 
becomes 



ay = [l-{l + aK)-{K+ -K-)IN^ 



part J 



(9) 



The value of ax is not one in isospin asymmetric colli- 
sions, however, we will use aK ~ 1 because the effect of 
a non-unitary value is reduced by {K'^ — K^)/Npart 
which is a small quantity for collisions we consider. It 
should be noted that only one unit of strangeness per 
multi-strange hyperon is counted in Eq. (^). This is safe 
because the production of multi-strange hyperons is rel- 
atively small even at the SPS energies p^ . 

The AGS E859/E866 proton and deuteron measure- 
ments cover a broad rapidity range exploiting the (near) 
symmetry of the Si+Al and Au+Au systems |23]. We 
take the {dN/dy)p weighted average of (/(y))jv"^ i^i 
Eq. (|l|) over this rapidity range to obtain (/)^'". The 
unmeasured d/p ratio in the mid-rapidity region is ap- 
proximated by the dotted lines in Fig. ^ connecting the 
last point (i.e., closest to mid-rapidity) with the refiected 
one. For Si+Au we assume that the unmeasured d/p 
ratio at the more forward rapidities is the same as the 
last data point. We assign a 8% systematic error on 
(Z)^'" for Si+Au due to this extrapolation by examin- 
ing a smoothly dropping d/p ratio from the low rapid- 
ity data points. From (/)]^"^, we extract (/)§'" using 
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Eqs. @, and @. For all three systems » 1 and 
ay ranges from 1.03 in peripheral to 1.13 in central colli- 
sions. For Si-f-Al and Au-I-Au, we assign a 5% systematic 
error on due to the assumptions made to extract 

(/)§'" 9,nd errors on the correction factors. For Si-|-Au, 
we assign a 10% systematic error on (/)§'". 

Figure I shows K+fR- § and (/)|'" at the AGS as 
a function of A'part- Both quantities are nearly constant. 
It should be noted that the constant values do not nec- 
essarily demonstrate that / K~ and {f)^™ are corre- 
lated. However, it will be shown that they are correlated 
in central collisions at different energies. Therefore, it is 
fair to argue that they are correlated also in the same 
collision system and are independent of the collision cen- 
trality. Because K^/K^ is fixed at chemical freeze-out, 
the results may imply that the values at kinetic 

freeze-out, which presumably happens later than chem- 
ical freeze-out in heavy-ion collisions, are connected to 
the chemical freeze-out values in a way that is in- 

dependent of centrality. 

Since the difference between the kaon yields {K^ —K^) 
enters into as a correction, one may be concerned 

about auto-correlation. However, such concern is un- 
necessary because (1) no correlation is present between 
/K~ and K'^ —K~ (which increases with A'pait more 
than linearly ||l[), and (2) the — K~ correction to 
is small (at most 13%). In fact, one may correlate 
jK^ with the average nucleon phase-space density 
which is shown as the dashed curve in Fig. ^. 

It should be noted that the correlation exists only 
between the rapidity integrated values of (/}§'" and 
K+/K-- locally in rapidity, (Df'^y) and K+/K-{y) 
are not or much less correlated (especially in central 
Au-|-Au collisions). This may be due to longitudinal flow 
which can be different for K~^, , and nucleon [p4| . 
Similar problems are encountered in equilibrium model 
studies of particle ratios: thermodynamic relations be- 
tween various particle ratios are only valid for ratios of 
integrated yields; differentially they do not hold because 
of longitudinal and transverse flow psf . 

Now we study the correlation between K^/K^ and 
(/)§'" i'^ central collisions at different beam energies. 
The values of K~^/K~ and (/)^'" extracted from ex- 
perimental data are tabulated in Table ||. The /K~' 
values for S-l-S and Ni-|-Ni |^ are at mid-rapidity. 
The values for S+S || and Pb-|-Pb ||,^ are 

determined from d/p at mid-rapidity. These values are 
approximated as for the total yield ratios with the sys- 
tematic errors quoted in Table | The K+ / R- and 
values for Au+Au [0Jl|] and Ni-fNi [|,|o) are obtained 
from measurements at slightly different beam energies. 
For Ni+Ni we have assumed p/n = 1. 

Figure]^ shows / K" as a function of (/)^'". Both 
quantities decrease with beam energy. The data from 
various collision systems at different beam energies follow 
a similar dependence. Motivated by equilibrium models. 



we fit the data to power-law function, yielding 

/ / f\Kin \ 1-20±0.10 

f^ = l+ ^Ub ] (10) 

K- V0-0076± 0.00097 

The fit result is plotted in Fig. || as the dashed curve. 

The KaoS Collaboration observed a K'^ /K~ ra- 
tio in Ni-|-Ni collisions at the same equivalent (sub- 
threshold) beam energy that is about seven times smaller 
than in p+p near threshold. The KaoS Collabora- 
tion contributed this difference to a drop in the 
mass in Ni+Ni collisions, and inferred a mass drop of 
27{)tfo MeV/c^ by using the measured kaon excitation 
function |^. Since is a measure of matter den- 

sity, we argue that studying / against (/)^'" may 
provide a different and probably more direct approach to 
search for medium effect on kaon mass: If the K~ mass is 
modified by medium at low energy but not (or differently) 
at high energy, then the low energy K'^ /K~ should de- 
viate from extrapolation of the high energy data. This 
Letter is a first attempt of this new approach: The fit 
result in Fig. ^ seems to imply no difference between the 
SIS and AGS/SPS data; fitting the data excluding the 
Ni+Ni point to power-law function, yielding the result 
shown in the dotted curve, seems to give the same impli- 
cation. 

Moreover, as noted, the average baryon phase-space 
densities at chemical and kinetic freeze-out may be differ- 
ent; the more suitable variable to correlate with / 
is the chemical freeze-out one. It has been found in 
Ref. ||l6|] , under the assumption of chemical and ther- 
mal equilibrium, that the chemical freeze-out is 
twice the kinetic freeze-out at both the AGS and 

SPS. Using the thermal model parameters obtained from 
the Ni+Ni data @, the {f)'^ value is found to be the 
same as the value. In Fig. | if the {f)^ values 

were plotted, then the Ni+Ni data point would lie above 
a sensible extrapolation from the AGS and SPS data by 
about a factor of two. This is even more contradictory to 
K~ mass drop, as mass drop would give a reduced 
K+/K- ratio. 

However, one should be careful to conclude from the 
above results that there is no K~ mass drop in the Ni+Ni 
data. This is because: 

1. Although the power law parameterization is moti- 
vated by equilibrium models, it is not a priori for 
the relation between /K~ and (/)§'". In fact, 
data from C+C collisions at the SIS do not fol- 
low the parameterization shown in Fig. |^: there a 
large /K~ ratio of 39 ± 6 was observed; no d/p 
data is available, but it can be conjectured that the 
(/)!'" value is lower in C+C than in Ni+Ni. It is 
possible that the C+C data follow a different sys- 
tematic consistent with the AGS/SPS data, while 
the Ni+Ni data lie below this systematic. This 
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would be an interesting result because it implies 
a K~ mass drop in Ni+Ni collisions. 

2. Since the leverage of the AGS and SPS data is 
small, the systematic errors on these data are crit- 
ically important to the extrapolation to the SIS 
data. 

In order to address these concerns, more data are needed 
to fill in the gap between AGS and SIS. These data are 
underway from the E895 Collaboration |3^. New data 
from the SPS at 40 AGeV would be also valuable as they 
fill in the gap between the full-energy AGS and SPS data. 

In summary, we have extracted the average baryon 
phase-space density at kinetic freeze-out from the 
deuteron to proton ratio in Si-|-Al, Si-|-Au and Au-|-Au 
collisions at the AGS 18|. The average baryon phase- 
space density at kinetic freeze-out is found to be nearly 
constant over the collision centrality, which offers an ex- 
planation for the puzzling centrality independence of the 
charged kaon total yield ratio [Q. The correlation be- 
tween the charged kaon ratio and the average baryon 
phase-space density is studied for central heavy-ion colli- 
sions of various systems over a wide range of beam energy. 
It is found that the charged kaon ratio increases strongly 
with the average baryon phase-space density, and both 
variables increase with decreasing beam energy. It is ar- 
gued that studying the correlation between these vari- 
ables may provide a different and more direct approach 
to search for medium effect on kaon mass. Although no 
conclusion can be drawn from the present study regard- 
ing medium modification to the kaon mass in Ni+Ni col- 
lisions at the SIS Q , new data from the AGS low energy 
collisions should shed light on this question. 
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FIG. 1. The spatial average nucleon phase-space density at 
kinetic freeze-out, (/(j/))jv'°i as a function of the proton cen- 
ter-of-mass rapidity {y) in central and peripheral heavy-ion 
collisions at the AGS. Collision centralities arc indicated 
by the cross section ranges. The data points at backward 
rapidities {y < 0) are measured, and are reflected about 
mid-rapidity {y = 0). Errors shown are statistical only. The 
experimental systematic effects largely cancel in the d/p ra- 
tio. The solid lines connecting the data points are to guide the 
eye. The dotted lines are an extrapolation to the unmeasured 
regions. 
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FIG. 3. The / K~ ratio as a function of the average 
baryon phase-space density at kinetic freeze-out in cen- 

tral heavy-ion collisions. Errors shown arc statistical and sys- 
tematic errors added in quadrature. The data are tabulated 
in Table I. The dashed curve is a fit to power-law function. 
The dotted curve is a similar fit but excluding the Ni-|-Ni data 
point. 
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FIG. 2. The K+ /K' ratio (filled symbols, left scale) and 
the average baryon phase-space density at kinetic freeze-out 
(/)!'" (open symbols, right scale) as a function of the to- 
tal number of participant nucleons (A/part) at the AGS. Er- 
rors shown are statistical for K'^ /K~ and statistical and sys- 
tematic errors added in quadrature for The experi- 
mental systematic effects largely cancel in the / K~ and 
d/p ratios. Note that both quantities are nearly constant 
over centrality. The dashed curve shows the average nucleon 
phase-space density at kinetic freeze-out (right scale) ; the er- 
rors are comparable to the (/)§'" ones. 
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TABLE I. The K+ / ratio and the average baryon phase-space density at kinetic freeze-out 
ratios) in central heavy-ion colhsions. The factors qj, ajv, and cty take into account, respectively, 
asymmetry, antinucleon and (anti)hyperon yields. For K'^ / K~ and the first error is statistical, 

(in percentage) are introduced by the extrapolation in rapidity, the uncertainties in a/, ajv, and ay, and 
of some of the data. The experimental systematic effects largely cancel in the / K~ and d/p ratios. 



" (extracted from d/p 
the effects of isospin 
The systematic errors 
the preliminary nature 



Collision Beam Energy Centrality 
System per Nucleon (mb) 



K+/K- 



Q-N 




s + s 

Pb-FPb 
Si-FAl 
Si-FAu 
Au-I-Au 
Ni-HNi 



200 GeV 
158 GeV 
14.6 GeV/c 
14.6 GeV/c 



160 
340 
100 
260 



11.1, 11.6° GeV/c 420, 470" 
1.8, 1.93" GeV 100 



1.58 ±0.09 ±5% m 

1.83 ±0.02 ±5% ||^ 

4.45 ±0.32 |i| 

5.25 ±0.42 |i| 

6.35 ±0.48 |i| 

21.1 ±3.4 ±5% |2| 



0.0039 ± 
0.0069 ± 
0.0221 ± 
0.0269 ± 
0.0297 ± 
0.0879 ± 



0.0007 ± 10% 
0.0010 ± 10% 
0.0017 ±5% 
0.0006 ± 10% 
0.0009 ± 5% 
0.0021 ± 5% 



I LI 
I LI 



° The first number is for the /K data; the second is for the d/p data. Quoted for Ni+Ni are the kinetic beam 
energies. The Pb+Pb ay value is used because the total yield — is not available for S+S. 
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